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TECHNICAL MEMORANDUM X-73364 


MSFC SOLAR SIMULATOR TEST PLANE 
UNIFORMITY MEASUREMENT 

SUMMARY 


The Marshall Space Flight Center (MSFC) solar simulator lamp array Is 
designed to produce a uniform Illumination on a 1.2 by 2.4 m (4 by 8 ft) test 
plane. Prior to completion, an analysis of the simulator was made to predict 
the test plane uniformity, A procedure was also developed to measure the 
actual uniformity produced by the simulator after it was in operation. That 
procedure is described here and the results of the first measurement presented. 
The measurement revealed the test plane illumination was not as uniform as the 
analysis predicted. The best actual uniformity obtained had a standard deviation 
of 4 percent and a peak-to-pealc variation of -ill percent. The computer program 
used to analyze the measurement data is discussed, and a listing of the program 
along with example runs is given. The reason for the variation in uniformity 
from that predicted is not fully resolved, but seems to be due to variations in 
individual lamp flux levels in tlie actual array from the levels when measured 
in a special test box. 


I. INTRODUCTION 


A large low cost solar simulator has been constructed at MSFC for the 
evaluation of flat plate solar collectors. The simulator consists of 405 tungsten- 
halogen projector lamps with Fresnel lenses arranged in a 15 by 27 array. The 
Fresnel lenses are used to provide a near collimated beam of illumination. It 
will be used to test solar collectors with surface areas up to 1. 2 by 2.4 m 
(4 by 8 ft) . The simulator is located in Building 4019 and v/111 be part of the 
Solar Test Facility under development at MSFC. 


An analysis of the uniformity of illumination expected in the test plane 
was described in Reference 1. Now that the simulator is in operation, the actual 
uniformity has been measured. This report gives the procedure used to meas- 
ure the test plane uniformity as well as the first measurement results. 


II. DISCUSSION OF TEST APPARATUS 


A. Test Philosophy 

The teat plane measurement apparatus was developed to produce data on 
the vmiformity in scan lines. This was done to ease comparison with the com- 
puted scan lines given by the computer program reported in NASA TM X-64991. 
The computer program was set up to calculate the irradiance at 5 cm {2 in, ) 
Intervals across the narrow axis of the test plane. Scan lines at any desired 
spacing along the long aj.is can be calculated. The mechanical scanner was 
designed to produce the same type scans, except the output is a continuous 
signal recorded on chart paper. Figure 1 is a layout of the measurement 
system. 



Figure 1. Layout of test plane uniformity 
measurement apparatus. 
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B. Scanner 

The scanner is a motorized platform which has a 10 cm (4 in* ) square 
plate on which a detector can be mounted. It is driven by a small 11!5 Vac 
400 Hz synchronized motor which is geared to move the platform across the 
width of the test plane in approximately 30 s. 

The scanner moves on a rail which can be positioned at 15 cm (G in. ) 
intervals along the frame. The frame is approximately 1.5 by 3.0m (5 by 
10 ft) so that an area larger than the actual test plane can be scanned. The 
scanner bar is positioned in the frame by means of pins placed in holes. The 
holes are drilled in the frame along the long axis, and near the end the holes 
are placed at 5 cm (2 In. ) intervals to allow higher density of scans near the 
edge of the test plane. 


C. Detector Circuit 

The detector used on the scanner platform for the first series of tests was 
a Motorola MRD 300 phototransistor. The circuit used with the phototransistor 
is shown in Figure 2. Tlie output is recorded on a Mosley model 680 strip 
chart recorder. The collector base junction of the phototransistor is used as 
a current source into the operational amplifier to provide a linear output voltage 
for varying levels of Incident radiation. The linearity of the detector was 
checked by comparing its output with readings from an Eppley blacii and white 
pycnometer model 8-48. The output is very linear with a constant voltage on 
the tungsten-halogen lamps and using neutral density filters to reduce the light 
level. However, when varying the voltage of the bulb to reduce the light level, 
a difference between the two detectors was noted. Figure 3 shows a plot of the 
phototransistor output compared with the pyranometer, as the lamp voltage is 
changed. For comparison, a curve is also sho\vn which gives the expected 
change if the light output varied as a function of the square of tiie voltage. 

These data are normalized to 100 Vac. The phototransistor agrees with the 
pyranometer for higher voltages but drops off at lower voltages. The reason 
for the drop off is that the pyranometer has a better infrared response than the 
phototransistor. At lower voltages the peak output of the lamp is shifted 
slightly toward the infrared. 
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D. Motor Drive Control 


The scanner motor control cii'cuit is shown in Figure 4 . The circuit 
automatically stops the scanner at the end of the scan. Microswitches on the 
scanner are used to indicate the end of the travel. When the microswitch hits 
the contact at the end of the bar, power to I’elay 1 as well as relay 2 is cut off, 
but the capacitor Cj on relay 2 holds that relay in momentarily. The result is 
that a brief reversing pulse is supplied to the motor to serve as a brake. Power 
to the motor is cut off by means of relay 2 before the motor can actually revei’se. 
To restart the scan in the other direction, the reversing switch is changed and 
the reset button is pushed. The reset button has to be held down until the micro- 
switches are clear of the contacts. The red Indicator bulb is used to show that 
the end of travel has been reached. 


III. COM PUTER AN ALYS I S PROG RAM 
A. Background 

The data recorded by the tost plane measurement equipment are a voltage 
level on chart paper representing the light intensity level for scans across the 
narrow axis of the test plane. Sample points have to be manually read from the 
chart paper to calculate the statistics and present the data in a useful form. A 
computer program, Statistical Analysis for the Solar Simulator (SASS), was 
developed to analyze the data. The sample points at regular intervals along 
each scan are read from the chart paper and recoi'dcd in a computer file for 
storage. The computer px’Ogram reads the data from the file and performs the 
data analysis. 

The number of points recorded from the scans and the total number of 
scans used can be varied. The normal procedure used was to read 29 points on 
each scan and use a total of 19 scans. Twenty-nine points on a scan at 5 cm 
(2 in.) intervals cover a range of -71 to 71 cm (-28 to 28 in.). Nineteen 
scans cover the test plane plus 15 cm (6 in.) on each end. Within the computer 
program the extra points and scans arc dropped to provide data only on the test 
plane. The test plane area can then be defined as the best 1. 2 by 2. 4 m (4 by 
8 ft) area in the recorded data. The program also allows any number of points 
or scans to be dropped in case the collector to bo tested is smaller than the 
standard test plane. 
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Figure 5 Is a simplified flow diagram of the compute, orogram. The 
full Fortran listing is given in Appendix A. The program is \ ritten is XDS 
extended Fortran IV for use on an XDS Sigma V time share computer. Except 
for the plot subroutine, the program can be run on any I'ortran computer with 
very little modification. 


6. Data Input 

Table 1 presents a listing of the input section of the computer program, 
lines 1 to 37. A listing of the variables used in the input and the main section 
of the program and their definitions is given in Table 2. 

The dimension statement on line 3 sets the size of the arrays. The 
arrays are now set for a limit of 27 scans with 31 points on each scan. By 
clianging the dimension statement, the limit can be changed to any value. 

These data are read in from a file on unit 100, When running the program, 
the input file name has to be assigned to unit 100, Appendix B gives a sample of 
an input file named T2, Appendix C gives an example run of the program 
together with the set commands to inform the computer which file will be read 
by unit 100. Line 9 is a rewind fiie statement which assures that the read 
statement will start with the first line in the input file. Line 11 is a read state~ 
ment to input the title comment. The format statement, line 16, allows up to 
80 characters in the title. Lines i? to 22 are the commands to input the control 
variables. Comment statements are included to define each variable. The 
foimat statements 110 and 111 allow a general input; the numbers can be fixed 
or floating point and only have to be separated by commas. The actual sample 
data points are Input by statements 23 to 34. The points are input into the 
B(J, K) array, where J Indicates the points on the scan and K indicates the 
scan number. 

These data are stored in the input file with consecutive lines containing 
the same point on each scan; for example, the first point on lines 7 to 37 of the 
input file T2 In Appendix B made up one scan. The second point on each line 
comprises the second scan, etc. Because of the limit on tlie size of a line in a 
file, all the scan points could not be put on one line; therefore, the scans were 
split up and NLPl number of them is read in first with statement 27. The rest 
is read in with statement 34, One line in a file will hold 72 characters; 
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TABLE 1. INPUT SECTION OF PROGRAM SASS 


li>uX ,§,iiCTIOJi OF FHCGfiAM SASS, 


■n 1-37 


- 

1.000 

C 

PKOGHAH SASS 

2.000 

C 

FrtOGhAK TO CALCJLATG STATISTICS GIJ DATA 

3.000 


DIl-iSriSION f; ( 3 1 , 27 ) , 0 ( ? 1 , 27 ) , IIATA ( 20 ) 

il.OOO 

C 

lOATA IS COMl'ib'ilT AhhAi 

5.000 


DlKtiliSIui^: X( 3 1 ) , ASi OS{ 3 1 ) 

6.000 

C 

i; 

7.000 

G 

IHPUT SECTION 

6.000 

C 


9.000 


REWIND 100 

10.000 

C 

READ Iri DATA 

11.000 


HEAD (100,100) IDATA 

12.000 


READ ( ICO, 1 10) ilSCAN , NFOS , i'PCOUF 

13.000 


HEAD ( 100, 111) SSL0C,TSi; C 

1^1.000 

c 

SSLOC-POSITION OF FIRST SC Ah 

15.000 

c 

TSIi'!C-Thc, SCAM IliChlM'iEiiT DISTANCE 

16.000 

100 

FORMAT (20A4) 

17.000 

c 

NSD-!iUMBER SCAN TO START TEST PLANE DATA 

16.000 

c 

NED-WUHPER SCAN TO HhO TEST PLANE DATA CM 

19.000 


KEAD( 100, 1 1 1 ) USD, NED 

20.000 

c 

MPDS-NUMBER POINT TC START DATA SCANS Oi. 

21.000 

c 

NPDE-NUMBER POINT TO END DATA SCANS ON 

22.000 


READ( ICO, 1 1 1 ) HPDS.NPDE 

23.000 


READ (100,110) NL,NLP1,AMV 

24.000 

110 

FORMAT (5G) 

25.000 

111 

FORMAT (20) 

26.000 


DO 200 J=1,NL 

27.000 


READ ( 100,250) X ( J ) , fILP 1 , ( fi ( J , K ) ,K= 1 , 1'LP 1 ) 

28.000 

200 

CONTINUE 

29.000 

250 

rORMAT (1G,NG) 

30.000 


IF(NLP1 .GE. NSC AN) CO TO ’DO 

31.000 


READ (100,111) ML,NLP2 

32.000 


MPsNLP1+NLP2 

33.000 


DO 300 J = 1,NL 

34.000 


READ ( 100,250) X ( J ) , I'LP2 , ( b ( J , ?0 , R=NLP 1 + 1 , NP ) 

35.000 

300 

CONTINUE 

36.000 

C 


37.000 

C 

END OF INPUT SECTION 


REFliOiHJCIBILlTY OF THfl 
PAGE IS POOBJ 
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TABLE 2. LISTING OF VARIABLES USED IN THE INPUT 
AND THE MAIN SECTION OF PROGRAM SASS 
(IN THE ORDER OF THEIR APPEARANCE) 


Variable 

Definition 

B(J,K) 

Input data array 

D{J,K) 

Array containing points of test plane 

IDATA{I) 

Comment array to store title 

X(I) 

Array containing the position of each point on 

- 

the scan 

ASPOS(I) 

Array containing the scan locations 

NSCAN 

Number of total scans 

NPOS 

Number of points on each scan 

NPCODE 

Code to print or plot input data 

SSLOC 

Location of first scan 

TSINC 

Distance between scans 

NSD 

Number of scan on which to start test plane 

NED 

Number of scan on which to end the test plane 

NPDS 

N;unber of point on scan to start test plane 

NPDE 

Number of point on scan to end test plane 

NL 

Number of lines for this input statement 

NLPl 

Number of points on line for this input 
statement 

AMV 

A multiplication value to normalize the data 

NLP2 

Number of points on line for this input 
statement 

NDUM 

Number of points on each scan used in test 
plane 

ND 

Number of scans used in test plane 

ADUMA 

Dummy variable to hold scan position 
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consequently, the number of digits used for a scan point will determine the 
maximum number of scans on a line. As in file T2, typically 11 scans are put 
in the first section. 


C. Main Section of Program 

Table 3 presents a listing of the main section of the computer program 
SASS. In statements 41 to 43 the data are normalized. The value of the vari- 
able AMV is selected to make the maximum value 100 percent, or it can be 
selected to convert the detector output to power units such as W/ cm^ or 
Btu/ h-ft^. The proper value for AMV to convert the output to power units is 
obtained by measuring the same point on the test plane with the phototransistor 
and the pyranometer. 

Lines 45 to 47 reduce the data array to the desired test plane and store 
the new points in the D(J,K) array. The size test plane would depend on the 
size collector under test. For evaluation of the simulator, the maximum test 
plane size of 1.2 by 2,4 m (4 by 8 ft) is used. Lines 53 to 54 define the vari- 
ables NDUM and ND, where NDUM is the number of scans in the test plane 
and ND is the number of points on each scan which lie in the test plane. State- 
ments 51 to 59 print out the statistical data on the test plane. All the actual 
computations take place in the subroutine ARMS which is described In Section 
m.D. 


The scan locations are set up in statements 60 to 63 and stored in the 
ASPOS(l) array. Statements 64 to 73 use NPCODE to select the print or plot 
options. If NPCODE is set equal to zero, no print or plot of the input sample 
data will take place. If the variable NPCODE is set equal to one, the input 
scaii data will be printed out along with a printout of the points used in the test 
plane evaluation. If NPCODE is set equal to two, the Input scan data \vill be 
plotted but no printout of the data will be obtained. No option is set to obtain 
both a printout and a plot as it is Impractical on the CRT display terminal. 
The printout is usually set to be produced on the line printer. Sections m. E 
and m. F describe the SPLOT and DPRINT subroutines and explain how to set 
the output device code to obtain the desired output. 
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table 3. MAIN SECTION OF PROGRAM SASS 


X'lAIN SECTION OF j>R0GHAM S«5 
TY 39-77 ^ 

39.000 C START INITALIZATION 

itO.OOO C 

111.000 DO 500 J = 1,.\,'P0S 

112.000 DO 500 K=1,NSCAN 

43.000 500 B(J,K) = J.fOJ^AMV 

44.000 C REDUCE ARRAY TO TEST PLANE SIZE 


45.000 

46.000 

47.000 1000 

48.000 C 

49.000 C 

50.000 C 

51.000 

52.000 2200 

53.000 

54.000 

55.000 

56.000 C 

57.000 

58.000 2400 

59.000 

60.000 
61.000 
62.000 

63.000 3000 

64.000 

65.000 
66.000 
67.000 

68.000 3100 
69.000 
70.000 
71.000 
72.000 

73.000 4000 

74.000 

75.000 

76.000 C 

77.000 C 


DO 1000 J = NPDSjKPDE 
DO 1000 K = NSD.K'ED 
D(J-(NPDS-1) ,K-(.iSD-l ) ) = b(J,K) 

PRINT STATISTICS ON ThST PLANE 

WRITE (108,2200) IDATA 

FOHMATC 1L1 ,20A4,//, 'TEST PLANE DATA*) 

NDUK = i^jpDE-RPDS+l 
ilD = NED-NSD+I 
CALL ARHS(D,iVD, NDUL) 

PRINT STATISTICS CM FULL nRASuREMSNT DATA 
WRITE(106,2400) 

F0HF;AT( 1H0, 'FULL FEASURLM^.NT DATA SET*) 

CALL ARMS { H , USCAN , MFCS ) 

ADUMA s SSLOC + TSINC 

DO 3000 I = 1,NSCAM 

ADUHA 5 ADUMA - TSIiJC 

ASPOS(l) s ADUMA 

IFOiPCODii.EC.O) CO TO 4000 

IF(!'!PCODE.N£. 1 ) CO TO 3100 

CALL DPRINT(X,B,ASP0S, 1 , NSCAN , 1 , NPOS , IDA'l A ) 

CALL DPhIMT ( X , b , ASPOS , USD , hED , liPDS , NPDF, IDATA ) 

IF(NPC0DS.NE.2) GO TO 4000 

CALL TPAUSE 

CALL SFLOT ( X , B , MPOS , hSCA'i , IDATA ) 

CALL TPLOT(0.0, 170.0,0,0) 

CALL ARMS(D,ND,NDUi'0 

CONTINUE 

STOP 

END 

« 

END OF MATH SECTION OF PROGRAM 


eepeoducebility of the 
PAGP i IS POOR 
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D. Subroutine ARMS 


Subroutine ARMS computes the statistical data. The subroutine calculates 
five values to characterize the uniformity. 

RIVK — Root mean square value computed from a zero intensity 
reference. 

FTP — The difference between the highest and the lowest sample point 
value. 

EBAR — The mean value of all the data points. 

SIGMA — The standard deviation of the measured values from the mean 
value. 

VIvIAX — The largest data point measured. 

A listing of subroutine ARMS is presented in Table 4. Table 5 presents 
a list of the variables used in the subroutine. The statistical values calculated 
by the program are slightly different in form from the standard approximation 
formulas. The standard approximation formulas as stated in Reference 2 are 
given by equations { l) , { 2) , and ( 3) ; 


EBAR = 


n 



n - 1 


RMS 



n - 1 


( 1 ) 


( 2 ) 


SIGMA = 



(y. - y)^ 


Va 


n - 1 


( 3 ) 
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TABLE 4. SUBROUTINE ARMS OF PROGRAM SASS 


tx 60-130 
80.000 
81.000 
62.000 
83 . COO 

84.000 

85.000 
86 . COO 
67.000 
88.000 

69.000 

90.000 

91.000 

92.000 
93.000 

94.000 

95.000 

96.000 

97.000 

96.000 

99.000 

100. coc 
101.000 
102. OCO 

105 . 00 c 

104.000 

105.000 
106.000 
1G7.00C 
lOd.COO 
1C9.000 
IIO.OOC 

111 . coo 
1 12.000 
1 13.000 
1 14. coo 
1 15.000 
1 1'j .000 
1 17.1)00 
1 IO. CoC 

1 19 . coo 

120. CC0 

121. COC 

1 22 . vOO 

125 . occ 

1 24. coc 

125 . COO 
1 2u . COC 

127.000 

126. COO 
129. OCC 
1 50 . COO 


SUBKOuTUvS ARMS(G'1,H,N) 

SUBROUTIMg AHllS IS TO COMPUTE STATISTICS 
or IS SAMPLE DATA ARRAY 
02 IS ARRAY OF SQUARED VALUES 
G3 IS ARRAY OP DEVIATION SQUARED VALUES 
N IS NUMBER OF POINTS IN EACH LINE 
M IS NUMBER OF LINES IN ARRAY 
DIMENSION G1(31,1),G2(31.27),G3(31,27) 

COMPUTE TAP, TOTAL NUMBER OF POINTS lU ARRAY 
TAP = FLOAT (ti»H) 

CALCULATE MEAN VALUE-SBAR 
CALL AFUNC(G1,EBAR,M,N) 

FIND G2 AND G3 ARRAYS 
DO 100 J = 1,N 
DO 100 K = 1,M 
G2(J,K) = G1(J,K)»«2 
100 G3(J,K) = (G1(J,K)-EBAR)*»2 

DU'M s 0.0 

CALL AFUNC(G2,DUM,M,U) 

RMS s SQRT(DUH) 

CALL AFUHC(C3,DUK,M,N) 

SIGMA s SORT(DUH) 

FIND MOST NEGATIVE NoMBER 
C = G1(l,1) 

DO 300 J s 1,W 
DO 300 K s 1,M 
IF(C-G1(J,K)) 300 , 300,290 
290 C = G1(J,K) 

300 CONTIMUE 

C IS NOW MOST NEGATIVl -DUMBER 
FIND HOST POSTIVE KUi -ER 
D = G1 ( 1 , 1) 

DO 400 J = 1,N 
DO 400 K = 1,M 
IF(G1(J,K)-D) 400,400,390 
390 D s G1(J,K) 

400 CONTINUE 

u IS NOW MOST POSTIVb ;UIM5BR 

PTP = D-C 

'.‘i.AX = D 

Fill NT GUT DATA 

1C00 FGKMATd HO, 'STATISTICS 0;i ARRAY DATA*) 

WhITEC 108 , 1100 ) 

11CG FORMAT (1 HO, 3X, ’ RMS * , 1 1 X, ’ FTP • , 1 1X, *EBAR*,9X, 
X'SIGHA' ,8X, 'V>iAX') 

VI rilTiri ( 1 C 8 , 1 200 ) RMS , FTP , EBAF , SIGMA , VMAX 
120C FORMAKIH ,2(F8.4,5X),F1C.5,2(5X,F8.4),/) 

10b, 1400) TAP 

1400 FORMATdHO, *NUHBEB OF POINTS IF ARRAY = *,2X,I4) 
RE'TuRfi 
EN D 
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TABLE 5. LIST OF VARIABLES IN SUBROUTINE ARMS 


1 

Variable 

Definition 

G1(I) 

Sample data array 

G2(I) 

Array of squared values 

G3(I) 

Array of deviation of squared values 

N 

Number of sample points on each scan 

M 

Number of scans 

TAP 

Total number of sample points 

DUM 

Dummy variable 

RMS 

Root mean squared value 

SIGMA 

Standard deviation 

C 

Smallest sample point value 

D 

Largest sample point value 

PTP 

Difference between D and C 

VMAX 

Redefinition of D 


where n is the number of sample points, is the value of ith sample point, and 
y is the mean value of all sample points (EBAR). 

It should be noted that slightly different equations are used when working 
with sample points which do not represent a function. The sample points are 
then considered a population, and n is used instead of n-1 (see Reference 3). 

The formulas given in equations (l), (2), and (3) are approximation 
formulas used for sample points representing a one-dimensional function. The 
irradiance of the test plane is a two-dimensional function where the irradiance 
is a function of both the x and y directions in the test plane ( Fig. 6) . To under- 
stand the equations that will be used to calculate the mean value (EBAR) , root 
mean squared value (RMS) and the standard deviation (SIGMA, it is instructive 
to review the derivations of equations (l), (2), and (3) . 


15 



Figure 6. Representation of irradiance 
function F(x,y). 


The proper definitions of EBAR, RMS, and SIGMA as given in Reference 
3 are given by equations (4) , (5), and (6) : 


EBAR = — ^ 
b - a 


b 

/ F(x) dx 
a 


RMS 


1 

b - a 


/ F**(x) dx 


V2 


SIGMA = 



V2 


f (F(x)-EBAR)^ 


( 4 ) 


(5) 


( 6 ) 



All three of the previously mentioned formulas can be written as one equation 
in the form of equation (7) : 


A = 


1 

b -a 


b 

/ G(x) dx 
a 


(7) 


In equations (4), (5), and (6) the function F(x) is integrated over the range 
from a to b, and in equation (7) the function G(x) is integrated over the range 
from a to b. The function G(x) is assigned the following values in order to 
give equations (4), (5), and (6): 


E G(x) = F(x) then A = EBAR 

EG(. , = F®(x) then A = (RMS)2 

E G(x) = (F(x) - EBAR)* then A = (SIGMA)* 


Since all three equations are similar, the derivation of equation (l), using 
G(x) = F(x), will at the same time give equations (2) and (3). 

E F(x) is a function (as represented in Fig. 7), then a set of sample 
points also can represent the function (as shown in Fig. 8), In terms of G(x). 
Performing the integral from a to b is the same as finding the area under the 
curve from a to b. When working with a set of sample points, the area is 
approximated by adding the area segments of width Ax under each sample point 
as shown in Figure 9. The formula for the integral is given liy equation (8) : 


n 

G(x) dx y y.Ax . (8) 

i=l " 


KJ 

/ 
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—♦I Ax — - 

Figure 8, Sample representation of function G(x). 

Since Ax is a constant, the summation can be rearranged as in equation (9) 


n n 

Z yM = Z • 





Figure 9. Approximation of area under function G(x) 
in region a to b. 


From Figure 9, it can be seen that Ax is given by equation (lO) ; 


Ax = 


b ~ a 
n - 1 


Substituting the value ot Ax into equation (9) gives equation (ll) : 


Ax 


i=l 


b - a 
n - 1 


n 



If this value for the integral is put Into equation (7) , then the value of A is 
given by equation ( 12) , 



A 




f 


( 12 ) 


1 b ~a 
b - a n “ 1 


n 

Z y, 

1=1 


which reduces to equation (13), 


A 


1 

n - 1 


n 


Z 

i=l 




(13) 


Equation (13) is the same as equation (l); therefore, we have derived 
equations (2) and (3) as well. For the limit, as n approaches infinity, the 
summation is truly equal to the integral. However, for small values of n, 
the approximation can be improved. By referring to Figure 9, it can be seen 
that half of the Ax wide area for the first and last point is outside the range 
of a to b. If one half of the area for each point is dropped, equation ( 8) is 
rewritten to equation (14) . That is, 


b n-1 . y Ax 

/ G(x) dx J y Ax + + ““ — 

a i=2 


(14) 


or rearranging to equation (15), 


/ g(x) dx Ax 


n-1 (yi+y„) 

2 * ~2 

i=2 


(15) 


The value of Ax is the same as in equation (lO) ; therefore, A is now equal to 
equation (16); 


A c- 


n-1 


n-1 (yi+y„) 

I n 

+ — ~ 


Z yi 

i=2 


( 16 ) 
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The improvement In the approximation of equation (16) over equation 
(13) is not very substantial for a one-dimensional function with a practical 
number of n sample points. However, the improvement can be important when 
worldng with a two-dimensional function. 

The two-dimensional ease is derived using the same technique as for the 
one-dimensional case. The two-dimensional formulas equivalent to equations 
(4), (5), and (6) are given by equations (17), (18), and (19); 


EBAR 


1 

(b - a)(d - c) 


b d 

f J F(x,y) dxdy 
a c 


RMS 


1 

(b - a)(d - c) 


b d 

/ / F^(x,y) dx dy 
a c 


(17) 


(18) 


SIGMA = 


1 

(b - a)(d - c) 


b d 

/ / (F(x,y) - EBAR)2 dx dy 
a c 


(19) 


where a and b are the limits of the function F(x,y) along the x axis and c and d 
are the limits along the y axis. Again, using a general function G(x,y) to 
derive all three values at one time, the double integral is approximated by a 
summation of AxAy volume sections centered on each sample point z,,^ (Fig. lO). 
The formula for A would be given by equation (20) ; 


jk 


A 


1 

(n - l)(m - 1) 


n m 
j=i J*' 


( 20 ) 


where n is the number of sample points along the x axis and mis the nuniber of 
sample points along the y axis. 


SAMPLE POINT 


• — Ax f— *J 

NOTE: DASHED LINES ENCLOSE THE ACTUAL AREA 
UNDER THE FUNCTION BY DEFINITION. 

Figure 10. Approximation of area under function G(x,y) . 


As in the one-dimensional case, the approximation can be improved for 
small numbers of n and m if the sample points fall on the edge of the function. 
If the volume outside the area bounded by a, b, c, and d Is discarded, the 
resulting approximation is given by equation (21) : 


b d (n-1 m-1 (z + z. ) 

/ / G(x,y)dxdy « AxAyj ^ ^ 

a c j=2 |li=2 


, m-1 + z, ) , fm-1 

1 V , 11 Im' , 1 V 
















The values of Ax and Ay are given by equations (22) and (23) : 


Ax = ^ (22) 

n “ 1 

Ay = . (23) 

m - 1 


Because the right-hand side of equation (21) is cumbersome, the following 
variables are defined: 



n-1 

m-1 

VI 

= E 

Z 


j=2 

k=2 


m-1 


V2 

= E 

k=2 

^Ik 


m-1 


V3 

II 

^nk 


(z.i + z ) 

jk 2 


(z-. + z- ) 

' 11 Im 


(z . + z ) 
nl nm 


(24) 


(25) 


(26) 


Using equations (21-26) , the value of A is given by equation (27) : 


(Vl + i V2 + I V3) 
“ (n-l)(m- l) “ 


(27) 


Ih the computer program the values of EBAR, RMS, and SIGMA are 
calculated by finding the three different values of the function G(x, y) . Three 
sample point arrays are used. The first, labeled Gl, contains all the sample 
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points as measured and is used to calculate EBAR. The second, labeled G2, 
contains all the sample points squared and is used to calculate the RMS value. 
The third, labeled G3, contains the squared values of the deviation of all the 
sample points from the mean value EBAR. A subroutine is called from the 
ARMS subroutine to calculate the value of A for each case. The subroutine is 
called AFUNC and is described in Section m. E, 

After calculating EBAR, RMS, and SIGMA, subroutine ARMS continues, 
calculating the lowest sample point, C, then the highest, D, and finally the 
difference which is labeled FTP. Lines 121 to 128 print out all the values 
along with the total number of points, TAP. 


E. Subroutine AFUNC 


A Fortran listing of the subroutine AFUNC is presented in Table 6. 

The variables used in the subroutine are given in Table 7. The subroutine is 
used to approximate the integral of a two-dimensional function represented 
by sample points where the outside sample points fall at the limits of the area. 
It actually computes the height of a box which if multiplied by the area under 
the sample points would be the integral of the function. 

Pour variables are transferred through subroutine AFUNC. The first, 
AVI, is a two-dimensional array containing the sample points. The second, 
AV2, is the value of A in equation (27) of Section ni.D. Tlie third, M, is the 
number of points in the y direction. The fourth, N, is tlie number of points 
in the x direction. 


F. Subroutine DPR I NT 

Subroutine DPRINT is used to print the sample point values. A listing 
of the subroutine is presented in Table 8. A list of the variables used is 
given in Table 9. The subroutine is called twice in the main program, once to 
print all the input sample points after normalization and again to print only the 
points used in the test plane. 
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SUBROUTINE 
'n 134-161 

134.000 

135.000 C 
13t.000 C 

137.000 C 

136.000 C 

139.000 

140. OOC C 

141.000 

142.000 

143.000 

144.000 
14S.0C0 

146.000 

147.000 

146.000 

145.000 

150.000 

151. COO 

152. COO C 

153.000 

154.000 

155.000 4C0 

156.000 

157.000 

156.000 C 

159.000 

160.000 
161.000 


TABLE 6. SUBROUTINE AFUNC OF PROGRAM SASS 

AFUiJC OF 


A 


100 

200 


300 


PltOGrtAK SASSj 

SUBROUTIi^E AFui.C(AV1,AV?,H,,<) 

SUBFOUl 1:.'L iC CALCcLriifc- liCTthOAL 
INTL’RGAL IS VULIKE Uiiijfch AKFA GiytM 51' AVI Ahuh 
RESULT AV2 IS HEIGHT OF VCLU:'-^' f CX 
AVI IS AFh'AY OF SAiiFLf 
DlijtNSICIi A VI (3 1,1) 

Calculate vi 
VI s 0 

uO 200 J = 2,i.-1 

100 K : 2 ,i .-1 

= vuAvi(j,i:^ 

= V1 + (AV1(J, 1 )+AV1(J,,;) )/?.o 


FOI.'ITS hfPEHSEi:Tlr;G G(X,Y) 


DO 

VI 

VI 


COt-iPUTL V2 
V2 r 0 

DO 300 fC = 2,. -1 
V2 r V2+AV1(1,K) 

* V|4 (AV1(1,1)+AV1(1. ;;))/?. r 
= .5*V2 
CCLPoTe V3 
V3 5 O.G 
D(J 400 IC r 2,F'-1 
V3 = V3+AV1(.i,K) 

V3 = V34(/;V1(i., 1)4AV1( ,,f.))/?.C 
V3 s 0.5*V3 

CALCULATE hEIOfil WHICH TIMES A}cA CIVES VOLU’’F 
AV2 = (VUV24V3)/FLCAT(( -l)5:(i .1)) 

F.ETUft'j 

END 


TABLE 7. LIST OF VARIABLES IN SUBROUTINE AFUNC 


Variable 

Deflnitioii ! 

AVl(J.K) 

AV2 

M 

N 

VI 

V2 

V3 

Sample point data array 
Result of Integral calculation 
Number of scans 
Number of points on the scans 
First term of integral 
Second teim of integral 
Third term of Integral 


reproducibility of ma 
19 unnii 
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TABLES. SUBROUTINE DPRINT OF PRCXSRAM SASS 


SUBRCUTlUt DPhIhT Ol-' SkSsY 

U 163-179 

1C 3 . COG SUbPOUTI :i£ DFHINT ( F , ft , .Sr , .-SS , G£S , , IJEF , 10 ) 

161*. COO GliltJiSIOu A(31,l),i5F(l),IL(l ),?(!) 

165. OCO WHlTtdOt, 100) (10(1), I = 1,20) 

166.0C0 100 Fun!iAT(lill,20AlO 


U7.000 
16b. CCO 
169. COO 
170. OCO 
171. OCO 
172.000 
173. OCO 
171 *. OCO 

175.000 

176.000 
177. COO 
176. COO 
179 . CCO 


/iPITr:(lC8,110) 

no FOniiAT(1r!O,»0ATA»,?1X,»SCA:’ FCGIT lOM • ) 
ilS = fJc.'£-r'S£+l 

Uhl'l 1:1 ( 1 C d , 12 c ) ii.S, (SF(I ) , I 2 \J-S,i'tS) 

120 rOnriATdH ,» f Ui IX,: (l- v . 1 , 1/ ) ) 
iihITf.dOfe, 150) 

150 FOr.GATdi! ,119(11-)) 

00 500 I = i;3F,r;nF 

snn£(1C6,510) F(I ) , l'55 , ( A ( I ,.J ) , J = •''^3,.jh.S) 
500 COirilOUh 

510 FOt:!:AT(H» *I * , IX, '.(F5. 1 , IX) ) 

HrX^U^rt 

SiJ:.- 



TABLE 9. LIST OF VARIABLES IN SUBROUTINE DPRINT 

I 



Variable 

OeRnltion 


P(I) 

Location of sample points on scan 


A(J,K) 

Array of sample data points 


SP(I) 

Location of scans 

^ _ Ktr . 

NSS 

Number of scan to start printing 


__NE8_— 

— Number of scan to e«d printing 

^ .... 

NSP 

Number of point on scan to start printing 


NEP 

Number of point on scan to end printing 


IA(I) 

Comment or title array 


NS 



The subroutine is almost standard Fortran, but it should be noted that 
a variable format statement is used to control the number of sample points 
printed. Some computers cannot handle a variable format statement. 

Using the set commands, the printout can be assigned to different devices. 
Usually, the unit 108 is set to the line printer if PCODE is assigned the value of 
1 for a printout. 


G. Subroutine SPLOT 

Subroutine SPLOT is used to plot the scan sample data in a three- 
dimensional format as an aid to visualizing the test plane uniformity. Table 10 
presents a listing of the subroutine and Table 11 gives a listing of the variables 
used. The plot package used is a special package developed by Computer 
Science Corporation for MSFC. It was developed for a Datacraft computer and 
made available on the Sigma V, Since this subroutine is computer-dependent, 
it will not be explained in detail. If the program is used on a computer other 
than the Sigma V, a new plot routine %vill have to be developed. 

The set command is used to control the location of tlie printed data in 
the plot routine. All write statements in subroutine SPLOT use unit 150 and 
the command set F:150 - ME is used when the NPCODE is assigned the value 
of 2 for a plot of the input data. 


IV. RESULTS OF FIRST MEASUREMENTS 


The computer program described in Reference 1 was used to calculate 
the uniformity of the 405 lamp solar simulator. Given the assumption that all 
lamps have the same intensity level at the same voltage, the program gave the 
results shown in Figure 11. The test plane showed only a tendency to decrease 
in Intensity near the edge. Using the procedure developed by Lewis Research 
Center, 405 lamps were tested and sorted for placement in the array in a 
manner to provide the most uniform test plane intensity. Using the measured 
values for the lamps in the computer program, the uniformity was again cal- 
culated and the results are shown in Figure 12 . The standard deviation 
increased from 1.09 to 1.45 percent. The PTP value was still less than 7 per- 
cent which was substantially less than the tolerance placed on the unifoimity, 

±7 percent deviation from the mean value. 


TABLE 10 . SUBROUTINE SPLOT OF PROGRAM SASS 


TY I 6 O -236 

180.000 


181.000 


182.000 


163. COO 


154.000 


185.000 


166. COO 


167.000 


166.000 


169.000 


190.000 


191. OOL 


192. COC 


193.000 


194 . COO 


195.000 

100 

196. COO 

lie 

197.000 


198.000 


199.000 


200. LOO 


201. OCO 


202. LOO 

16 G 

2IJ3.0C0 

150 

204. OCO 


205. CCO 


206. COG 


207 . CCC 


206. OCO 


209 . OCO 

200 

2 10. OCO 

2ic 

211.000 


212.000 


> 13. COO 

220 

214.000 



SUBROUTINE SPLOKTX/rPI.riX.ijJf, IDATA) 

DIMEtSSION TX(31),TPI(31,27),ICATA(20),1A(10) 

DATA (IA(I),I = 1,10)/1HP,1he,lHR,1HC,lH£,1HK,1H7 
XlriA, IhG, 1HE/ 

CAlL BEGIN (300) 

CALL ERASE 

CALL SCALE (6. 0,3. 0,200.0, 150.0) 

CALL VECTOR 

CALL AXIS (0.0, 0.0, 60.0, 100.0,6.0,20.0,1,1) 

CALL ALPHA 

HA = -36;X = -10.0;Y = -10.0 
DO ICC I = 1,11 
KA = MA+ 6 ;X = X+6.0 
CALL TPLCT(X,Y,0,0) 

v:rite( 150 , no) na 

COi'iTlilUE 

cOnMAT(I3) 

X = -18.0;Y = 100.0 
DO 150 I = 1,10 

1 = Y-8.0 

CALL T?L0T(X,Y,0,0) 
v,F.ITE(150, 160) IA(I) 

FOS:,AT( 1A1 ) 

COKTIflUE 

X = -12.0;Y = -21.0;NA = -20 
DO 200 I = 1,6 
i-;A = NA+20;Y = Y+20.0 
CALL TFLOT(X,Y,0,0) 

NHITE(150,210) NA 

COiJTIfiUE 

FOH^■AT(Iif) 

CALL 7rLQT(0. 0,-20. 0,0,0) 

WRITE( 150,220) 

FOr.i-iATC DISTANCE PROM CEi'ITERLIWE - INCHES*) 

CALL VECTOR 


» 


215. COO 

2 16. C CO 
217.0l' 0 
21&.0GC 
2i5.rrc, 

2 2 C . u 0 0 

221 . U.0 

222. LUO 

2 2 3. LOO 
22 A. LUO 
225. LOG 
226. OCO 
227. LUG 
226.000 
229.000 
230.000 

231 . OtO 

232 . C CO 

233 . OLO 
PS'l-OLL 
236. COO 

236 . CCC 


ZERO = 30.Ci;XSHIPT = 3.0;Y3HIPT = 5.0 

DEL = rX(l)+ZERO 

XS = -XSfiIFT;YS = -YSKIFT 

DO 500 I = 1,wY 

YS = YS+Y SHIFT 

XS1 = (XS+Ds.L+XSHIFT)«(1.0-FLOAT(I-1)*.C1) 
CALL TFL0T(XS1,YS,0,0) 

XS = XS+XSHIFT 
DO 250 J = 1,iiX 

X = (TX(J)+ZER0+XS)«(1.0-FLOAT(I-1)».L1) 

1 = 'iVI{J,r)+Y3 
CALL ThLOT(y, Y, 1 , 0 ) 

25c COLTIl uE 

CALL TPL0T(X,Y3, 1,0) 

300 COMTIMLE 

CALL TPLOT(C.C, 160.0,0,0) 

CALL ALPHA 

'■vklTE( ISC.'lOO) (IDATA(I),I = 1,20) 

JiCO F0H:-IAT(2CAli) 

CALL TPALSE 

hETuRil 

ti.'D 
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TABLE 11. LIST OF VARIABLES IN StJBROUTINE SPLOT 


/ 


Variable 

Definition 

TX(I) 

Array containing location of sample points on scan 

TPI(J,K) 

Data array of sample points 

NX 

Number of sample points on scan 

NY 

Number of scans 

IDATA 

Comment or title array 

IA(I) 

Axis label array 

NA 

Variable used to number the axis 

X 

Variable used to position number of the x axis 

Y 

Variable used to position number of the y axis 

ZERO 

Shift necessary to position data in center of plot 

XSHIFT 

Shift in X direction for each scan for three- 
dimensional effect 

YSHIFT 

Shift in y direction for each scan for three- 
dimensional effect 

DEL 

Shift of sample point to position correctly on 
X axis 

XS 

Temporary variable for shift value in x direction 

YS 

Temporary variable for shift value in y direction 
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COMPUTED TEST PUNE UNIFORMITY - ALL UMP8THE SAME - 14 ft 


RMS PTP EBAR SIOMA 



UMAX 

100.0000 


Figure 11, Plot of test plane uniformity with all lamps 
set at the same value. 

TEST PLANE UNIFORMITY WITH MEASURED LAMPS 



DISTANCE FROM CENTERLINE (In.) 


Figure 12. Plot of computed test plane uniformity 
vitli measured lamps. 



On May 20, 1976 the simulator was turned on and the actual test plane 
uniformity measured. A plot of the results is shown in Figure 13. After 
reviewing the data, it appeared the test plane measurement apparatus had not 
been centered under the simulator; therefore, a second test was conducted. 

The second measurement in ( Fig. 14) gave results that are now considered 
typical. The FTP value was 22. 6 percent (±11. 3 percent), and the standard 
deviation was 3. 6 percent. 

Measurements mth the test plane at greater distances than the standard 
2.74 m (9 ft) distance from the simulator show that the edge drop-off increases 
as the computer program predicted. Figure 15 is a computer plot of the test 
plane at 3, 66 m (12 ft) , and Figure 16 is a plot of the test plane as measured 
at the same distance. Figures 17 and 18 give similar comparison for 4. 57 m 
(15 ft). 


At this time it is not clear why the measured results show a much larger 
variation in uniformity than that predicted by computer analysis. At first, it 
appeared that light spillover from one lamp to the adjacent Fresnel lens might 
be causing the nonuniformity, but a test was conducted in which an aluminum 
honeycomb 2. 3 cm (0. 9 in, ) thick with holes approximately 0.5 cm (0.2 in. ) 
in diameter was placed between the lamps and the Fresnel lens. The intensity 
was greatly reduced with the honeycomb (approximately 80 percent), but the 
uniformity was not significantly changed when normalized. (The honeycomb 
was painted black to reduce reflection in the channels. ) A check of the individ- 
ual lamp intensities were made while in the full array, and they were not as 
measured prior to placement in the array. They were placed correctly by 
serial number, but the intensities were different. Why they were different is 
unexplained. It could be due to such factors as contact resistance in the lamp 
sockets or voltage drops in the power lines feeding the individual lamps. 

A subsequent test of lamps in the special test box has revealed that dif- 
fei«nt intensity levels, on the order of 15 to 20 percent, will be measured if a 
lamp is rotated 180 degrees in the lamp socket. Tlie reason for this is unex- 
plained also, but . could account for some of the variations of the lamp intensities 
when measured in the actual array. 

After operating the simulator. It was noted that the Fresnel lens sagged 
due to the heat, resulting in a warped surface with a curve up to 6 mm (0.25 in. ) 
deep across the 12, 7 cm (5 in. ) diameter, perpendicular to the two mounting 
screws at the edge. The warping seems to have little effect on the optical prop- 
erty, but could account for some of the test plane variation. The warping is a 
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TEST plane UNIFQtiMiTY with MEASURED UMPS 12 ft 


RMS PTP EBAR 5I0MA UMAX 

97.1384 13.8177 08.85237 3.6812 100.0800 


20 I h I I 

.30 -24 -18 -12 -8 0 6 12 18 24 30 

DISTANCE FROM CENTERLINE (in.) 



1 


i 


Figure 15* Plot of coniputed test plane uniformity | 

at 3.66 m (12 ft). ] 

, ^ 

SSTP TEST PLANE AT 3.6S m (12 ftj 

BMS PTP EBAR SIOMA UMAX ! 

B0JS03 zg,2375 90.22125 4,9314 999999 | 



Figure 16. Plot of measured test plane uniformity 
at 3.66 m (12 ft). 





problem in that thermal cycling may result, in physical failure of the lens. 

Dr. McCrlckard of ET44 has an effort underway to investigate the warping 
effect and tc develop means to minimize the problem. 

Another anomaly observed was that the flux level from the lamp array 
was about 30 percent higher than expected. TMs had to be corrected because 
the maximum life of the lamps is achieved only if the lamps are operated 
between 90 and 110 Vac. To I'educe the test plane intensity to Air Mass 2 and 
still operate above 90 Vac, screen wire was placed between the lamps and the 
Fresnel lens. The screen wire reduced the intensity by approximately 30 per- 
cent and left the imiformity unchanged. 


V. CONCLUSIONS 


The 405 lamp solar simulator has been placed into operation. The test 
plane uniformity is not as good as predicted but is within acceptable limits. 

The reason for the difference between the measured results and the predicted 
performance continues to be evaluated. It is felt at this time that the difference 
Is the sum of several effects, but that the primary reason is the variation in 
flux levels for the individual lamps wlien measured in the test box from the 
levels measured when placed in the actual simulator. Other possibilities 
include the warping of the Fresnel lens or individual manufacturing variations 
in the lens. 

The total flux level at the test plane was higher than the predicted value. 
Compensation is obtained by placing screen wire between the lamps and the 
Fresnel lens which reduced the test plane intensity level to the design value. 

Although the process to manually read die sample points from the chart 
paper is tedious, the procedure used to measure the test plane uniformity 
worked very well. 

The computer program to reduce the measurement data uses slightly 
modified formulas to calculate the statistics on the test plane uniformity. The 
formulas were modified to allow the edge of the test plane to be defined by 
sample points. The modifications are not required if the sample points are 
taken at the center of each sample area. 
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The test plane measurements at different distances from the simulator 
indicate that the computer program given in Reference 1 to calculate the teat 
plane uniformity gives accurate results; however, it will not be fully proven 
until a comparison can be made in which the individual lamp intensities are 
accurately known. The large variation observed in lamp performance in the 
simulator indicates that it is extremely difficult to control all the parameters 
to insure that the individual lamp intensities measured in a test fixture will be 
repeatable in the simulator operating conditions. 

To more fully understand the variations from the predicted performance 
and to provide operational flexibilily and efficiency, further evaluation and 
possible design improvements are recommended as follows: 

a. Test a variety of lamps to determine expected variation due to qualify 
control, voltage, orientation, operating temperature, etc. 

b. Build and test lamps in a large simulator in the lab to determine 
problems when going from the test of one lamp to an assembly of lamps. 

c. Conduct test on a large sample of Fresnel lenses to determine manu- 
facturing variations and changed in optical performance due to warping. 

d. Design and build an automated scanner to measure the test plane 
uniformity. An inexpensive microprocessor with a motor driven scanner would 
quickly provide the test plane uniformity or the illiunination on a collector of 
varying sizes, at a large saving in manpower and time involved to obtain the 
measurement results. 
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C PRBQRAM 8ASS 

~c pworaR' t6 CaVCOlate statistics ONI Data 

DIHENSION B<3l«87)4DOKS7),tDATA(S0) 

"c n?Atr i's comment' a'RRav 

DIMENSION X01),ASP0801) 

“C 1 

C INPUT SECTION 

“£ ♦ 

REWIND 100 

~C ■READ'TNTDAtA 

READ (iOO«tOO) tOATA 

RfeAD(ioo,iio> NscAN,"WBi;i^perbiyE“ 

REAOdOOiUl) SSL0C#TSINC 
"C BSE:TCiF8SlTlON OF PlRSf^CSFl 

c tsinc-the scan increment distance 

~To3 P5RHAT <20A4) 

C NSDuNUMBeR scan to start TEST PLANE DATA 

-C NED.NUMBER scan to end ' TeST“P1:XNe 'WtX'W>(~ 

L NPDS*NU °* ^R^OT^'^^^STARr-IJA- r SCTA1M5-SN 

c npoeInumbIr point to end data scans on 

READdOOdin "NPUSTNPOE 

READ (lOOtllO) NUjNLPIiAMV 

— TTO FORMA nWI 

111 FORMATOO) 

DO SOO Ji ITNC 

READ (100^250) X(J>,N{.P1/ (B( J#K)#K*1jN(,PI) 
—550^ CONTINUr 

2S0 FORMAT «18#NQ) 

rFt N LP l «a E*NSCANl ob ~to~^oo 

READ (1004 111) N1.4NLP2 

NPfNtPT+NtP2 

DO 300 J « IjNL 

RgAt)-n:o o4 g so ) "x ( - j n An;:prj-rr(OTn7iaiciW7Ni^ 

300 CONTINUE 

"C 9 

C END OF INPUT SECTION 

"C * 

C START IN IJAL I Z AT_I^ 

DO 500 Jil^NPOS 

00^00“ KWl/NSCAl^ 

500 B{J«K) f B(JjK)*AMV 

-c REDUCr ARRAY- T0"TBST^CTNE“3Ur* 

DO 1000 J 9 NPDSjNPDE 

OO-IOOO 'i?"l TJSDTNEO" " ’ 

1000 D (J* ( N PD8 *1)<K *(NSD »1 )) 9 B(JiK ) 

^ ^ STICS ON TES T PLaNe 

WRITE (10842200) IDATA 

~220^r TORMA TTIHI 4'?0 AJFi >/rrTe8T“PCANn5'ATAr i 
NDUM • NPDE«NPDSA1 


reproducibility oe the 
ORIGINAL PAGE IS POOR 
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NO « NeD-NSD«t 

CAL L - AW St07l<0>-mtWl 

C HPRlNT STATISTICS SN PUIU Mea8URe«ENT OATA 

wRrTrtio8*-2Aoon 

8400 FORMAT (IH0< 'FULL MEASUREMENT DATA 8et») 

C ALL A RM S<B#MSCANiN P8^ 

ADUMA « SSLOC * TSINC 

—OS -3000‘*T~wt> N8CAN 

ADUMA i ADUMA • TSINC 

- 3 0 00 "As pest n * a p u m a 

IRNPCeOEtEOfO) 09 TO 4000 

tffNPes-PE tW E il ) 09 TO' 3 10 0 

CALL DpRINT ( X iB# ASpOSi i#NSCAN« li NpOS« lOATA ) 

c all D PR I N T ( X/ B , a SP 9S<N SP < N£ B , N PO S/Np O m By 

3100 IF(NPCBOE«NEfS) 09 T9 4000 

C ALL TP AUS E 

CALL SPLST(X#BiNP9S»NSCAN4l0ATA) 


CALL ARMS(D«ND4NDUM) 

- 4 000 -C9 NTI -N0E 

STOP 

iNO 

c 4^^ or main oe c ti o n o r p noo nA M 

C WBR9UTINE ARMS 

subroutine ARMS(Q1«M«N) 

-c s ub ro u t ine A RMS-tS TO CO M PUTE" ST A TTST ICS 

C Qi IS SAMPLE DATA ARRAY 

-e Sr'l-S-AlWAY-OP-SOtJAftED^VAtyeS- 

C 03IS ARRAY OF DEVIATION S QUARED VALUfS 

C M IS number OF LINES IN ARRAY 

DiMENSFON etF3lil->*0gi3lia y)*Q3 < 3l<gy > 

C COMPUTE TAPj total NUMBeR OF POINTS IN ARRAY 

TAP f Ft O ATtN*M-)- 

C CALCUL ATE MEAN VALUE.EBAR 

C find os and (33 ARRAYS 

DO 100~3~f“l# N 

DO 100 K > liH 

os’tOTKr-i "tmoiXT w2 

100 Q3(JiK) • (Q1 (JjK)«EBAR)*«2 

DWM-"*~OVt) 

CAL L AFU NC(Og<DUM,M,N) 

CALL AFUNC(03#DUM«M«N) 

S IQ M A g-S Q RT fOU H-l 

C FIND MOST negative NUMBER 

C— f-Q-H 1 #-H 

DO 300 J f 1«N 

DO — 30 0 K B- ly M 

IF(C«01(J«K}) 300<300«290 

— S50-C-4 01id4 K > - 

300 CONTINUE 
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C C IS New MOST NEQATIVE number 
C PIND M6ST PSSTIVC NUMBER 
D 9 Qia«i) 

DS 400 J 9 1/N 
0» 400 K 9 1/M 
IP(G1(J/K)«0) 400i400«390 
390 D 9 Q1(J«K> 

400 continue 

C D IS New MeST POSTIVe number 
PTP 9 O-C 
VMAX 1 0 

c print eUT Data 

lOOO P0RMAT(lH0/ 'STATISTICS 9N ARRAY OAfAM 
WRITE IXOBi 1100) „ 

1100 PORMATnNo^3Xi f RMS i j Hxi f PTP ( / UXi iEBAR»# 9X^ 
X'SIGMA'jgX/ 'VMAXM 

WHITE <108/ 1200) RMS/FTP/EBAR^siOMA, VMAX 
1200 FORMATuH /2(F8.4/5xJ*P10*5#2(5x<P8#4)#/) 
WRITE<108/1400) TaP 

1400 P9RMAT(lHo/ 1 number 8F POINTS IN ARRAY • »/2X/l4) 
RETURN 
END 

* 

SUBROUTINE APUNC 
SUBROUTINE APUNCf AV1/AV2/M,N) 

subroutine to calculate INTERQAl 
INTeRQAL is volume under area QIVEN by AVI ARRAY 
Result avs is height of volume sox 
AVI IS array Of Sample points RepeRSentinq a<X/Y) 

DIMENSION AVK31/1) 

C calculate VI 
VI • 0 

DO 200 J 9 2/N*l 
DO 100 K e 2/M-l 
100 VI 1 V1*AVI(J/K) 

200 VI * VI*(aVI<J/1)*AVKJ/M))/2iO 
C COMPUTE V2 
V2 « 0 

DO 300 K B 2/M-l 
300 V2 i V24AVK1/K) 

V2 ■ V2+ (AVlU/l)*AVltl/M) )/2»0 
V2 t #5#V2 
C COMPUTE V3 
V3 8 OfO 

DO 400 K » 2/M.l 
400 V3 9 V3*AVt<N/K) 

V3 # V3*tAVl(N/i)*AVl <N/M) )/2.0 
V3 9 0«5«V3 

c calculate height which TIMES AR£A QIVeS VOLUME 

AVa • <V14V2+V3)/FL0AT< (N-1)#<M.1) ) 

RETURN 

END 


C SUBReUTlNg OPRIKT T8 PRlNj INPUT DATA 

SU1R9UT 1 NE OPR INT ( P# A# SP# NSS, NgS , N8P# NEP# I D ) 
DiMENSieN A(31il)»SP(l)4lO(l)iPm 
WRITEU08«100I (IDdlit 9 l#SO) 

100 P9RMAT(lHl480A4> 

wRTTErioajUoi 

110 FORMATUHO# »0ATA»j31X4 *SCaN POSITIRR’) 

NS • NES*NSS4>1 

WRITEdOSilSO) NSj(SP(n#l ! NSSiNES) 

ISO FORMAT fiw #♦ POS •flXjNTFSvmxn 
MR1TE{108«150> 

190 FORMATIIH dlRUH-M 
Oe 500 1 i nsp#nep 

NRITEdOdjSlOl Pn)4NS4(AU#J)iJ • NSS4NES) 

SOO CONTINUE 

fflO FORMATUH 4F4#0i »njlX4N(F5*liIXn 
RETURN 

eno 

subroutine SPUOTITX^TPI^NX/NV, IOATA) 

OIMENSKN TX(3ll4TPI(31<27WrOATA<ao)jIAdO» 

Data (IAU),! f l,lO)/lHP4lHE4lHR4iMC4lHEjlHN,lHT, 
XIHAi1HQ,iHE/ 

CALL BEG IN (300) 

CALL ERASE 

CALL SCALE<6tO/3tO«300«Oil90fO) 

CALL VECTOR 

CALL AXISlOfO/OtOiSOtOilOOtOjStO^SOtOjUi) 

CALL alpha 

NA 9 "SSiX 9 •10*0|Y 9 alOtO 

DO too I i i«n 

NA 9 NAA6JX 9 X«6i0 
CALL TPLOTfX*Y#04C) 

WRITEdSOillO) NA 
100 CONTINUE 
110 F0RMAT(I3) 

X 9 »18iO|Y • 100-*'0 
00 150 ! 9 IdO 
Y # Y»8.0 

CALL TPLOT(X#Y4040) 

WRITEdSOdOO) lAd) 

160 FORMATdAl) 

JTTO CONTINUE 

X 9 9lSiOiY 9 •SltOiNA 9 920 
00 200 I 9 1/6 
NA I NA«20iY 9 Y920tO 
CALL TPLOT(X/Y/0/0> 

WRITE|150j810) NA 
200 continue 
210 FORMATdA) 

CALL TPL0T(0t0/920«0/0/0) 

WRITEd50#280) 

220 FORMAT (» DISTANCE FROM CENTERLTNE • 2•NCHCSi^ 

CALL VECTOR 

EERO 9 30.OJX8HIPT v 3*0|Y9M1PT i S«0 
DEL 9 TXd)*2ER0 


X8 f iixSHIFTiVS • *VSHIFT 
00 300 I i 1»NY 
YS t YS+YSHIFT 

XSl • {XS+0£L+XSMrPT)*{l«0^FL9ATn.l)*,0l) 
CALL TPteT(XSljYS# 0 , 0 ) 

XS i XS*XSHIFT 
00 25o J • 1#NX 

X • (7X(J)+ZERa4XS»#(lt0-FL0AT(I-n*,0l) 

Y n TPKJjD+YS 
CALL TPL6T(X/Y,1,0) 

250 CSNTINUF 

CALL TPL0T(X,YS«1,O» 

300 continue 

call TPL8T(Oi0/l8OtO#O#O) 

CALL ALPHA 

WRITE{150/400) (I0ATA(I)#I • 1^20) 

400 format (20A4) 

CALL TPAUSE 

RETURN 

END 


44 


APPENDIX B 


SAMPLE OF INPUT FILE 
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E 


uF lut'Uf HL^ FCj> t'KOUHAi’ iAob FAGL 1_ 

Ti 1-l*0 

1. C00 cjSTF-G/ 2;V76-;.I;J1:- FX,2..fi i.u;:- 0 IC^TIO: .': 

2. G00 19,51,1 
3.00C tO,C 
n.oco 2, It 
5.000 4,2G 


6 . 

COC 

31 , 11 , 1.005 









• f 7 ^> t 


7 . 

OOC 

- 30 , 57 . 5 , 75 , 7 b. 5 , bO 

.”, 9.1 

f 92 

. 7 ,' 

li 

t 

65 . 

c- ,(. 

r 

7 t. 

t. 

ooc 

- 2 t , 56 , 72 . 9 ,b 1 ,c.i 4,0 

3 . 3,6 



.. i: 

■ •* f 

(. 

1.7 

,7r 

f 

Jl.n 


5 . 

COC 

— ^ 0,5 6 , 7 ^, ('£..(•, '.'b. 1 


, r. 

, f. 6 , 

• 


, i .2 

,■ ( 

• 


1.9 

10 . 

OLC 

— £; , 5 . (• , 7 ^ . '• , b< . 0 , (i 7 

• 3,99 

r 

• - » 

l 4 |. 1 

,96 

• 

^ r 

. 1 

f 

" 3 • If'*' 

11 . 

OCC 

- 22 , 5 l, 72 . 7 ,*“ 1 .^,b 7 




9 .( 

f 



7 

, r i 7 


12 . 

CCO 

— 2 C, 5 '~’.b, 75 .^,' ' .^-t 

‘ 

• , 

■ 

1 . * 

f 

j 

• • 

1 ,' 

** 

.( ,c: 7 , 

67.' 

15 . 

(.CO 

- 1 b, 59 . 5 ,T^),bC. 2 ,tO 

• t C 

,63 

. 7 , 4 - 

- • ' 

f 

c 4 . 

1 ,- 

X 

• f * 1 • 

1 , ' 


OCC 

- 1 u , 60 , 75 . 1 r bO . 3 , 0 2 

. 5 , <• (' 

. 6 , 

' ^ * 
v> ^ • 

, 9 ’* 

• 

4 ,4 

4.1 

f 

. 2 , >.• 

* f * 

15 . 

COC 

- 1 ^ 4 , 61 , 75 . >., 79 . ‘ 4 , r 2 

. 9,9 2 

, 6 o 

, 64 . 

1 , 

K 

r 

• - f 

6 w, 

< 

':.6,r. 


It. 

CCO 

-l?,(; 1 . 1 , 7 t', 7 ';,“ 2 .t 

, C' 3 . 5 

,65 

r - 

• .> 1 > 

6 . 6 

f 

: 6 . 

•• ,' 

7 

,>• 9 . 4 , 

'7 

17 . 

ccc 

- 10 , 02 . 6 , 77 . 0 , 79 . 7 , 

o 5 ,c 5 

,67 

. 7 ,< 

6.7 

f 

<• 1 . 

9,6 

0 • 

.4,91. 

7 ft *4 

1 b. 

000 

- 6 , 63 . 7 , 77 . 0 , 1 . 0 .b , 1 . 

v. 0,9 

c . 9 

f ^ ' • 

1,9 

h 

r 

• I 


7 

, 9 . J , 

\ 1 . 

19 . 

000 

— 6 ,oC.t.’, 4 ' 0 . l,cl. 9 ,>-' 

7 .c.,b 

1 ." 

,•^ 9 , 


c. 


,of. 

• 

5 ,^* 1 . - 

f ‘ 

20 . 

COc 

-A,fc. 7 .o,b 2 ,o 3 . 7 ,‘'v. 

7 ,'l. 

^'.> 

G > •* 

• ^ f 

r 


C 

t 

f 

4 - 1 . 7 , ' 

•7 

f 

21 * 

000 

- 2 , 66 . 0 , 63 , 05 . 7 , 90 . 

1 , 91 . 

9,9 

4 . 7 , 




• f 

r, 

3 . 9 ,br 

.7, 

22 • 

coo 

0 , o 9 . 2 , o 4 (,b 7 . 5 , 90,9 

1 • f L 

0 . ‘ 

r ^ 

f ' • 

f 

' I • 

r, 

'j 

. 1 , 

4*1 . 

r. 

■ ‘.'i ' 

, ^ • ■• , 

‘ - 

23 . 

too 

2 , 69 . c, 60. 1 ,(' 6 , jc,<. 


-6 , 

r».'> 

.?, 

C 

9,4 

'• f 




2 i*. 

ooc 

0 ,c 9 . 2 ,i 92 .< ,c•^ . 1 ,• 7 

• ^ f ' ’ ^ 

. 6 , 

•v 4 . c 

• r* 

I • 

f 

” 7 , 

•' 7 • 

r 

,'■‘6.6, 


25 . 

cot 

6,0 O ,b^ .i^,b0.t o.> 

f f ^ 

5 ,c 

5 . 6 , 

9 6 . 

r 

,64 

I 

• 

5 , ». . 6 


26 . 

ooc 

0 -, (• 7 . 7 , 01 . 5 , 95 . 5, 66 

• j , w w 

.6, 

b6,b 

0 . 9 

f 

bS 


1 

, ' ? • 6 , 

.1* . 

27 . 

,CGO 

10 , 67 . 3 ,bC , o 5 , 6tj . 3 , 

(-' 6 . ‘.i t 

o 7 . 

. 5,0 7 

• » f 

1 

b.b 


• 

6 , 94 ,/ 

7 • ' 

26 . 

, 0 t ;0 

12 ,o 7 . 7 , 79 .-", 6 't. 5 ,f 

f • 

6 ,t 

i'. . (•■ , 

•**. « ■ « 


, 9 '" 

c 

• ^ 9 

c 

V*' 

. . ^ 

^ ■ • 




. I , f <■ . 2 


30.000 I6,tc,t2.t .rT.T.iXt.’.f 1, '7. ?,f C .5 

31. OC Li 1 0 , 6 b . 2 , t f . 1 , > A , 1 , {■ J< , I , i.' , . , 1 , , I • , '■ . I. » f 0 , ^ 

32.000 2C , bt . '4 , c 2 . 0 , c . 2 , 0*^ . 1 , GO . 7 , ^ 7 , ‘ ' . ' , F ? , ^ ^ 0 . i) , r 7 

*•<3.000 22 , o7 , i) 1 . ? , o3 , >55 , Ot , 7 , b . 2 , C.2 , oO ,r( 

3 A . CCC 2U , b5 , 00 , 1 2 , c5 . ? , b t , c 6 . 3 , . 5 . 1 . 3 , 7' . f , i ^ 5 . '-• 

35.000 2o,bl.8,7G.t',79.*,^2 .i,f>3.7|b5.5»‘',-..'»'b.3»7*’,‘-2,' 


. 000 2o , 60 , 7^ , 7>.. . A , cO . t' , < < , ( ^ . 2 , ^ ,7 . 1 ,75.7 , f i' . ? , f 1 


' . 000 30 , 55. b , 70 . 5 , 72 , 7C , 79 , bO . 3 , 79 , 75 , 7^4 , 7b , 7^ . 2 
• .000 31 ,b 

; . GCC -30 , b1 . 5 , «;2 , 7H , 7G , Vo . ' , 74 , 7 1 . t , 5V . 1 
,'.0C0 — 2c,b2.^,< <'.^,bo,f O,< 1, / '•.o,7'-. 1»* 0.5 


kiKi*RODUClBILITY OF THfl 
ORIGINAL PAGE IS POOR 


w 


SflrIPLb ca-' Il^PUt flLii fOh t'l.uC.i’-V iAorf PftCk ?_ 

Ti iJ1-69 

it 1.000 -£6,L3,C3.i 1,a?,79.9,72.6,<''C.6 

i*2,000 -2il,63.i;‘,b5,iJ3,c3,a3,t0.i,7?,6l 
itS.OCC -22,6ft.6,fiG.b,t2.c,e5.9|f 3. 1^75,60.0 

44.000 -20,&4.7,t!5.7,i:.2.2,t4,7,c3.i',t1.f',77.4,60.6 

45.000 -r8,t2.4,04.6,{.1.4,c4,b2. 1,f C,77,62 

40.000 -16,l3.6,b4.?,c.1.b,b4.r,r'4.1,f0.5,7f:,5,CP.t 

47.000 -14,83,e5.7,b2.5,b5.3,b4.5,bC.2,77.7,f)3.2 

46.000 -12,&2,6,b6,C2.5,tD.2,a3.7,60,7i .6,f 4 

49.000 -I0,64,b6,62,4,c4.5,62.5,t 0. 1,76.7,^i^ 

50.000 -&,65,b7.6»62,64.*; ,t:2.7,'-0.2,7i .P,{ 3.'= 

5 1 . 000 -6 , 65 . 1 , eb . 4 , r,2 . 4 , 65 , 6 3 . t' , 6 1 . ? , 76 . 7 , 6 3 . 7 

52 . 000 -4 , 6 5 , 4, 66 , 62 . C , 64 . 7 , < 5 , i 2 . 5 , 79 , 0 ? . 4 

53.000 -2,66.5,67. 1 ,‘ 5. 2, <'5. 9, <-6.4, J.. 4, {-r, 03 

54.000 0, dC. 6, fco. 6, 04. 4, 67. 6, 67. 3,65.5,61.6,62 

55.000 2, 66. 5, 66, &4. 6, 67,7,67. 4, {'4, 62. 3, 02.3 
5 C . 000 4,67,60,65.5, 67. 7,05. 6 , 6 1 , :5 , 60 . 6 , (. 3 

57.000 6, 67. 2, 67. 4, 66.6,67. 2, 63. C,o1. 2, 79. 3, C2 

56.000 6,06.2,66.6,66.0,67.4,63,79.9,77,62 

59.000 10, 86. 4, 66. 9, 66, 66.6,63,76, 76, 61 

60.000 12,66.5,60. t, 65.6, t9, 4, 42. 6, 77. 2, 76. 1,59 

61.000 1 4, 66, 66. 0,64. 4, 66. 3, 61. 3, 76. 2, 73. 9, 5< 

62.000 16, 06. 1,66. 3, 64, 07.0,79.9,79, 74,96.0 

63.000 IS, 65, 64, 62, 3, 65. 5, 74,76.9,74.3,54.1 

64.000 2C,64.1,t-3,in.6,64.i ,76, 19.3,73.0,54.1 

65.000 22, 63. 2, 60, 6C, 63. 7, 79. 5, 7'''. '',73.9,54 

66.000 24,62.6,76.7,79.0,63.4,60.7,79.6,74.2,54.5 

67.000 26,63,76.6,76,62.6,40,79,75.4,95.2 

66.000 28, 62.5, 70. 4, 78. 3, 63.2, 60 . 4,tG. 7, 75, 55 

69.000 30, 79. 4, 74,7b, 80, 70.4,76. 5, 72, 54 

^ . . . . 


BBPEODUCBILH^ '^1' 
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DESCRIPTION OF INPUT FILE 


LINE 1 Coimnent or ntle input, (IDATA) can ocmtain up to 80 characters 

LINE 2 NSCAN, NPOS and NPCODE are Input on this line. 

NSCAN — Number of scans 

NPOS » Number of points on each scan 

NPCODE Code to print or plot input data. 

0 for no print or plot 

1 for printout on line printer 

2 for plot oa CRT 

UNE 3 SSLOC and TSINC are input on this line. 

SSLOC — The position of the first scan in inches, 

TSINC — The interval between scans in inches. 

LINE 4 NSD and NED are input on this line. 

NSD — Number of scan to start the test plane, 

NED — Number of scan to end the test plane. 

LINE S NPDS and NPDE are input on this line* 

NPDS — Number of sample point on the scan to start the test 
plane. 

NPDE — Number of sample point on the scan to end the test 
plane. 

LINE 6 NL, NLPl, and AMV are input on this line. 

NL — Number of lines in the first set of data. 

NLP1»-* Number of sample points on each line. 

AMV — Scale factor to multiply against each sample point. 

The line numbers may vary from here on depending on the number 
of sample points on each scan. 

LINES 7-37 These lines contain the first group cf sample points. The first 
lines start at the negative end and progress to the positive end. 
Each line contains the sample point on each scan at the same 
distance from the center point. 
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UNE 38 


KL and NLP2 are input on this line. 

KL Number of lines in the second group of data. 

NLP2 — Number of points on each line in the second group. 

LINES 39-69 These lines contain the rest of the sample data points. 
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APPENDIX C 

SAMPLE RUN OF PROGRAM 


-^B pmiK G SAGE Ntffl 
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EXAMPLE OF TO GET OUTPuT ON THt Llt^E PRliCTER^ 


CONTROL CGHiiAi:DS_ 


!SET F:1C0/T2 


JSET F:10G s LP 


I START TD2 
»STOP* 0 




t2fS7IO 


9Qi3Bl93 


3«*3£<» 100*0034 


•ttoiii 


41*4939 


97 181949 


4*810t 100*0034 




SSTP»6/a2/76*NI^E FTi 2M3 Rl^N- i<ie m(JDIFICATI»NS 


DATA 

P8S 

5<mo 

48*0 

42*0 

36*0 

SCAN P8SITI8N 
3p*o 24*0 18*0 

12»0 

6*0 

*0 

-6*0 

”12*0 

*18*0 

*24*0 

i 

til 

c 

• 

c 

»36«u 

•42*0 


77>5 

88*0 

93 *C 
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EXAMPLE OF RUN TO OBTAIN PLOT 
CONTROL COMMANDS: 

(SET F:100/ T2 (T2 is name of Input file. ) 

ISETF;150 = ME 

(EDIT T2 (Need to edit T2 to change plot code* ) 

*TY2 

2.000 19,31,1 
+IN2 

2.00019,31,2 

♦END 

(START TD2 (Command to run program SASS in file TD2. ) 

The following data Is output on tlie CRT 

SSTP-6/ 22/ 76-NINE FT, 2ND RUN- NO MODIFICATIONS 
TEST PLANE DATA 

RMS PTP EBAR SIGMA VMAX 

90.4249 22.5780 90,35193 3.6324 100,000 

NUMBER OF POINTS IN ARRAY = 425 

FULL MEASUREMENT DATA SET 

RMS PTP EBAR SIGMA VMAX 

88.0832 42.4935 87.81969 6.8102 100.000 

NUMBER OF POINTS IN ARRAY = 589 
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